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ABSTRACT 

The fluctuations in the emission of redshifted 21cm photons from neutral inter- 
galactic hydrogen will provide an unprecedented probe of the reionization era. Con- 
ventional wisdom assumes that this 21cm signal disappears as soon as reionization 
is complete, when little atomic hydrogen is left through most of the volume of the 
inter-galactic medium (IGM). However observations of damped Lya absorbers indi- 
cate that the fraction of hydrogen in its neutral form is significant by mass at all 
redshifts. Here we use a physically-motivated model to show that residual neutral 
gas, confined to dense regions in the IGM with a high recombination rate, will gener- 
ate a significant post-rcionization 21cm signal. We show that the power-spectrum of 
fluctuations in this signal will be detectable by the first generation of low-frequency 
observatories at a signal-to-noise that is comparable to that achievable in observations 
of the reionization era. The statistics of 21cm fluctuations will therefore probe not only 
the pre-reionization IGM, but rather the entire process of HII region overlap, as well 
as the appearance of the diffuse ionized IGM. 

Key words: cosmology: diffuse radiation, large scale structure, theory - galaxies: 
high redshift, inter-galactic medium 



1 INTRODUCTION 

The reionization of cosmic hydrogen by the first stars and 
quasars (e.g. Barkana & Loeb 2001), was an important mile- 
stone in the history of the Universe. The recent discovery of 
distant quasars has allowed detailed absorption studies of 
the state of the high redshift intergalactic medium (IGM) 
at a time when the universe was less than a billion years old 
(Fan et al. 2006; White et al. 2003). Several studies have 
used the evolution of the ionizing background inferred from 
these spectra to argue that the reionization of cosmic hy- 
drogen was completed just beyond z ~ 6 (Fan et al. 2006; 
Gnedin & Fan 2006; White et al. 2003). However, other au- 
thors have claimed that the evidence for this rapid change 
becomes significantly weaker for a different choice of density 
distribution in the IGM (Becker et al. 2007). Different ar- 
guments in favour of a rapidly evolving IGM at z > 6 are 
based on the properties of the putative HII regions inferred 
around the highest redshift quasars (Wyithe & Loeb 2004; 
Wyithe, Loeb, & Carilli 2005; Mesinger & Haiman 2005). 
However, Bolton & Haehnelt (2007a) and Lidz et al. (2007) 
have demonstrated that the interpretation of the spectral 
features is uncertain and that the observed spectra could 
either be produced by an HII region, or by a classical prox- 
imity zone. One reason for the ambiguity in interpreting 



these absorption spectra is that Lya absorption can only be 
used to probe neutral fractions that are smaller than 10~ 3 
owing to the large cross-section of the Lya resonance. 

A better probe of the process of reionization is pro- 
vided by redshifted 21cm observations. Reionization starts 
with ionized (HII) regions around galaxies, which later grow 
to surround groups of galaxies. The process of reionization 
is completed when these HII regions overlap (defining the 
so-called overlap epoch) and fill-up most of the volume be- 
tween galaxies. The conventional wisdom is that the 21cm 
signal disappears after the overlap epoch, because there is 
little neutral hydrogen left through most of intergalactic 
space. However, the following simple estimate can be used to 
demonstrate that 21cm emission should be significant even 
after reionization is completed. Damped Lya systems are 
believed to contain the majority of the neutral gas at high 
redshifts. Indeed, observations of damped Lya systems out 
to a redshift of z ~ 4 show the cosmological density pa- 
rameter of HI to be Ohi ~ 10~ 3 (Prochaska et al 2005). 
In the standard cosmological model the density parameter 
of baryons is Qt ~ 0.04, so that the mass-averaged neu- 
tral hydrogen fraction at z ~ 4 (long after the end of the 
HII overlap epoch) is f m ~ 0.03. This neutral gas does not 
contribute significantly to the effective Lya optical depth, 
which is sensitive to the volume averaged neutral fraction 
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(with a value that is orders of magnitude lower). However 
the redshifted 21cm emission is sensitive to the total (mass- 
weighted) optical depth of this neutral gas. Observations of 
the redshifted 21cm signal would therefore detect the total 
neutral hydrogen content in a volume of IGM dictated by 
the observatory beam and frequency band-pass. 

Although the 21cm emission after HII overlap is dom- 
inated by dense clumps of gas rather than by diffuse gas 
in the IGM as is the case before reionization is complete, 
we do not expect 21cm self absorption to impact the level of 
21cm emission. This conclusion is based on 21cm absorption 
studies towards damped Lya systems at a range of redshifts 
between z ~ and z ~ 3.4, which show optical depths to 
absorption of the back-ground quasar flux with values less 
than a few percent (Kanekar & Chengalur 2003; Curran et 
al. 2007). The small optical depth for self absorption is also 
supported by theoretical calculations of the 21cm optical 
depth of neutral gas in high redshift mini-halos (Furlanetto 
& Loeb 2002). Moreover, damped Lya systems have a spin 
temperature that is large relative to the temperature of the 
cosmic microwave background radiation, and will therefore 
have a level of emission that is independent of the kinetic 
gas temperature (e.g. Kanekar & Chengalur 2003). 

At 2 ~ 4, the brightness temperature contrast of red- 
shifted 21cm emission will be AT ~ 0.5mK. Moreover on 
the R ~ 10 co-moving Mpc (hereafter cMpc) scales relevant 
for upcoming 21cm experiments, the root-mean-square am- 
plitude of density fluctuations at z ~ 4 is a ~ 0.2. Hence, 
we expect fluctuations in the radiation field of ~ O.lmK on 
10 cMpc scales. The fluctuations in the 21cm emission signal 
after reionization are therefore expected to be only an or- 
der of magnitude or so smaller than the largest fluctuations 
predicted at any time during the entire reionization era (e.g. 
Wyithe & Morales 2007). In addition the sky temperature, 
which provides the limiting factor in the system noise at the 
low frequencies relevant for 21cm studies, is proportional to 
(1 + z) 2S , and so is a factor of ~ 3.4[(1 + z)/5] 2 ' 6 smaller 
at low-redshifts than for observations at z ~ 7. Thus, de- 
tectability of fluctuations in 21cm emission may not decline 
substantially following the overlap epoch, allowing the study 
of the entire HII overlap process in redshifted 21cm emission. 

In this paper we calculate the mean level of 21cm emis- 
sion through the overlap epoch to lower redshifts, as well 
as the statistics of the fluctuations about this mean due 
to patchy reionization and fluctuations in the density field. 
Throughout the paper we adopt the set of cosmological pa- 
rameters determined by WMAP3 (Spergel et al. 2007) for a 
flat ACDM universe. 



2 SEMI-ANALYTIC MODEL FOR 
REIONIZATION 

Miralda-Escude et al. (2000) presented a model which allows 
the calculation of an effective recombination rate in an in- 
homogeneous universe by assuming a maximum overdensity 
(A c ) penetrated by ionizing photons within HII regions. The 
model assumes that reionization progresses rapidly through 
islands of lower density prior to the overlap of individual cos- 
mological ionized regions. Following the overlap epoch, the 
remaining regions of high density are gradually ionized. It is 
therefore hypothesized that at any time, regions with gas be- 



low some critical overdensity Ai = Pi/(p) are highly ionized 
while regions of higher density are not. In what follows, we 
draw primarily from their prescription and refer the reader 
to the original paper for a detailed discussion of its motiva- 
tions and assumptions. Wyithe & Loeb (2003) employed this 
prescription within a semi-analytic model of reionization. 
This model was extended by Srbinovsky & Wyithe (2006) 
and by Wyithe, Bolton & Haehnelt (2007). We refer the 
reader to those papers and only summarise the basics of 
the model below. In the current work we limit our atten- 
tion to reionization due to population-II stars (Gnedin & 
Fan 2006; Srbinovsky & Wyithe 2006), which govern the 
final stages of reionization even in the presence of an ear- 
lier partial or full reionization by population III stars (e.g. 
Wyithe & Loeb 2003). 

Within the model of Miralda-Escude et al. (2000) we 
describe the post-overlap evolution of the IGM by comput- 
ing the evolution of the fraction of mass in regions with 
overdensity below Ai, 



Fm ( Ai) = f 'dAP v (A)A, (1) 
Jo 



where Py(A) is the volume weighted probability distribution 
for A. In a region of large scale overdensity 5 at z b s , the 
mass fraction Jm(Ai) [or equivalently Ai] therefore evolves 
according to the equation 



rfF M (Ai) _ 1 dn-fjS) R ( A n (i | $ D ( z ) \ dt 
dz no dz a :i c \ D(z b B ) J dz' 

(2) 

where D is the growth factor of linear density fluctuations, 
qb is the case B recombination coefficient, n c is the comov- 
ing electron density, -R(Ai) is the effective clumping factor 
of the IGM (see below), and dn-y/dz is the co- moving emis- 
sivity of ionizing photons. This equation was described in 
Wyithe & Loeb (2003), but is generalized here and below to 
regions of large scale overdensity 8 that differ from the aver- 
age IGM (Wyithe, Bolton & Haehnelt 2007). Integration of 
equation ((2)1 requires knowledge of TV (A). Miralda-Escude 
et al. (2000) quote a fitting function which provides a good 
fit to the volume weighted probability distribution for the 
baryon density in cosmological hydrodynamical simulations. 
This probability distribution remains a reasonable descrip- 
tion at high redshift when confronted with a more modern 
cosmology and updated simulations, although the addition 
of an analytical approximation for the high density tail of 
the distribution remains necessary as a best guess at cor- 
recting for numerical resolution (Bolton & Haehnelt 2007b) . 

Equation @ provides a good description of the evolu- 
tion of the ionization fraction following the overlap epoch 
of individual ionized bubbles, because the ionization fronts 
are exposed to the mean ionizing radiation field. However 
prior to the overlap epoch, the prescription is inadequate, 
due to the large fluctuations in the intensity of the ionizing 
radiation. A more accurate model to describe the evolution 
of the ionized volume prior to the overlap epoch was sug- 
gested by Miralda-Escude et al. (2000). In our notation the 
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appropriate equation is 

dQi 1 dn~f (S) 

dz n°_F M (A cr it) dz 



al. 1991) as 



a B (l + z) 3 R{A clit )n c [l + S 



-D(zobs) / dz 



dt 



dF M (A clit ) 



dz 



F M (A C 



(3) 



In this expression, Qi is redefined to be the volume filling 
factor within which all matter at densities below A c has 
been ionized, and R is the effective clumping of the IGM 

Within this formalism, the epoch of overlap is precisely 
defined as the time when Qi reaches unity. However, we have 
only a single equation to describe the evolution of two in- 
dependent quantities Qi and Fm (or equivalently A c ). The 
relative growth of these depends on the luminosity func- 
tion and spatial distribution of the sources. The appropri- 
ate value of A c is set by the mean separation of the ion- 
izing sources. We assume A c to be constant with redshift 
before the overlap epoch and compute results for models 
with values of A c = 20 and A c = 5. Prior to reionization, 
these values imply a mean-free path that is comparable to 
the separation between galaxies (Wyithe & Loeb 2003). Our 
approach is to compute a reionization history given a par- 
ticular value of A c , combined with assumed values for the 
efficiency of star-formation and the fraction of ionizing pho- 
tons that escape from galaxies. With this history in place 
we then compute the evolution of the background radiation 
field due to these same sources. After the overlap epoch, 
ionizing photons will experience attenuation due to resid- 
ual overdense pockets of HI gas. We use the description of 
Miralda-Escude et al. (2000) to estimate the ionizing photon 
mean-free-path, and subsequently derive the attenuation of 
ionizing photons. We then compute the flux at the Lyman- 
limit in the IGM due to sources immediate to each epoch, 
in addition to redshifted contributions from earlier epochs. 

We next describe our model for the emissivity of 
population-II stars. Following Wyithe & Loeb (2003) we as- 
sume the spectral energy distribution (SED) of population-II 
star forming galaxies, using the model presented in Leitherer 
et al. (1999). From this SED we calculate the number of ion- 
izing photons produced per baryon (-/V 7 ). Note that we ig- 
nore helium reionization in our modeling. We further assume 
that ionizing photons are produced primarily in starbursts, 
with lifetimes much shorter than the Hubble time, so that 
we may express the star formation rate per unit volume as 



dM 
~d~V 



r dF(S,z) 

/*— 77 Pb, 



(4) 



where pb is the co-moving baryonic mass density, and F(S, z) 
is the density dependent collapsed fraction of mass in halos 
above a critical mass at z. The factor /* (star formation ef- 
ficiency) describes the fraction of collapsed baryonic matter 
that participates in star formation. 

In a region of co-moving radius R and mean overdensity 
S(z) — S[D(z)/D(z b B )] [specified at redshift z instead of 
the usual z = 0], the relevant collapsed fraction is obtained 
from the extended Press-Schechter (1974) model (Bond et 



F(8, R, z) = erfc I 



5 C - 8(z) 



/2 ([a gal ] 2 - [a(R)] 



(5) 



where erfc(x) is the error function, a(R) is the variance of 
the density field smoothed on a scale R, and cr ga i is the 
variance of the density field smoothed on a scale R sa i- The 
latter scale corresponds to the minimum halo mass for star- 
formation, and both variances are evaluated at redshift z 
rather than at z = 0. In this expression, the critical linear 
overdensity for the collapse of a spherical top-hat density 
perturbation is S c ~ 1.69. 

Prior to the redshift of the overlap epoch, the relevant 
scale R in the calculation of the biased collapsed fraction 
(equation [5} is equal to the size of the region being con- 
sidered. However, following the overlap epoch the ionizing 
photon mean-free-path increases rapidly and quickly exceeds 
this scale. Once the mean-free-path is large, the ionizing pho- 
tons contributing to the background within a region of size 
R were not produced locally. At times when the mean-free- 
path A is greater than R we therefore assume that a fraction 
1 — 7? 3 /A 3 of the ionizing back-ground within a region of 
radius R was produced from sources in the mean-universe 
(we expect the contribution of the ionizing background from 
outside R to be only weakly correlated with the density of 
the region within R when A > The remaining fraction 
(J? 3 /A 3 ) is computed assuming biased star formation using 
equation (0 • Note that this dilution applies only to the ion- 
izing background, and not to the density of baryons within 
the region of size R. 

In a cold neutral IGM beyond the redshift of reioniza- 
tion, the collapsed fraction should be computed for halos of 
sufficient mass to initiate star formation. The critical virial 
temperature is set by the temperature (Tn ~ 10 4 K) above 
which efficient atomic hydrogen cooling promotes star for- 
mation. Following the reionization of a region, the Jeans 
mass in the heated IGM limits accretion to halos above 
Ti ~ 10 5 K (Efstathiou 1992; Thoul & Weinberg 1996; Dijk- 
stra et al. 2004) . We may therefore write the time derivative 
of the collapsed fraction 



dF 

dt veal 



(*) = 



dz 



dz 



dt v 



(6) 



where Q m is the ionized mass fraction in the universe. 

To describe the ionizing flux from stars we require one 
additional parameter. Only a fraction of ionizing photons 
produced by stars may enter the IGM. Therefore an addi- 
tional factor of / csc (the escape fraction) must be included 
when computing the emissivity of galaxies. Ciardi & Fer- 
rara (2005) include a review of existing constraints on / esc 
which suggests that its value is < 15%. More recent simu- 
lations (Gnedin, Kravtsov, & Chen 2007) show the average 
escape fraction from galaxies to be only a few percent over a 
wide redshift range. The product of the star formation effi- 
ciency and the escape fraction may be combined into a single 
free parameter (/*, eS c = /* x / OS c) to describe the contribu- 
tion of stars in our model. The ionizing photon emissivity 
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Figure 1. Models for the reionization of the IGM and the subsequent post-overlap evolution of the ionizing radiation field. In each panel 
two cases are shown, corresponding to different values for the critical overdensity prior to the overlap epoch (A c = 5, thin dark lines; 
and A c = 20, thick grey lines). We show the cases for the mean IGM with 5 = 0. Upper Left Panel: The ionization rate as a function 
of redshift. The observational points are from Bolton et al. (2007b). Upper Right Panel: The volume (lower curves) and mass (upper 
curves) averaged fractions of neutral gas in the universe. Also shown (dotted lines) is the fraction of the IGM yet to overlap (1 — Qi). The 
observational points for the volume averaged neutral fraction are from Bolton et al. (2007b), while the observed mass- fractions are from 
the damped Lya measurements of Prochaska et al. (2005). Lower Left Panel: The mean-free-path for ionizing photons computed using 
the formalism in §[2] The data points are based on Storrie-Lombardi et al. (1994). Lower Right Panel: The evolution of the mean 21cm 
brightness temperature (in mK) with redshift (solid lines). For comparison, the fraction of IGM yet to overlap (1 — Qi) is overplotted. 



may then be written as 

dn~, r , T dF . . 

^ = /W V 7 -(*)p bl (7) 

where pb is the co-moving mass-density of baryons. 
3 RESULTS 

Figure Q] shows example models for the reionization of the 
IGM and the subsequent post-overlap evolution of the ion- 
izing radiation field. Two cases are shown, corresponding 
to two different values for the critical overdensity prior to 
the overlap epoch (A c = 5 and A c = 20), with values of 
/*,csc = 0.0048 and / tjCSC = 0.0038 respectively. These cases 
do not represent the best fit to the data, but rather bracket 
the range of the overlap epoch redshifts (6 < z < 7) for 
which our model is consistent with observations at z < 6 
(without invoking an additional population of more massive 
stars at high redshift; e.g. Wyithe & Loeb 2003). The ex- 
amples have different values of A c , and thus illustrate the 
mild dependence of our results on this unknown parame- 
ter. We note that the values of /*, OS c required for our model 
to reproduce existing observations are in excellent agree- 
ment with external considerations. In particular our value 
of /*, C sc ~ a few xlO" 3 corresponds to product of recent 



estimates for the escape fraction (a few xl0~~ 2 ; Gnedin, 
Kravtsov, & Chen 2007), with estimates of the average star- 
formation rate (~ 10 _1 from the ratio between the mass den- 
sity in stars and baryons; Fukugita, Hogan & Peebles 1998). 

In the top left panel of Figure [T] we show the evolution 
of the ionization rate. The observational points are from 
the simulations of Bolton et al. (2007b; based on the obser- 
vations of Fan et al. 2006). In the top- right panel we plot 
the corresponding volume (lower curves) and mass (upper 
curves) averaged fractions of neutral gas in the universe. 
Also shown (dotted lines) is the fraction of the IGM yet 
to overlap (1 — Qi). The redshift where these curves drop 
to zero is normally quoted as the redshift of reionization, 
and these curves correspond approximately to the standard 
semi-analytic calculation (e.g. Haiman & Loeb 1997). How- 
ever our formalism allows for the calculation of both mass 
and volume averaged neutral fractions to lower redshifts. 
The observational points for the volume averaged neutral 
fraction are from Bolton et al. (2007b), while the observed 
mass-fractions are from the damped Lycv measurements of 
Prochaska et al. (2005) , and therefore represent lower limits 
on the total HI content of the IGM. Both curves show excel- 
lent agreement with these quantities, despite their differing 
by 3 orders of magnitude. In the lower left panel we plot the 
evolution of the ionizing photon mean-free-path. The data 
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e (arc-minute) 8 (arc-minute) 

Figure 2. Examples of auto-correlation and skewness functions. Upper Panels: The auto-correlation is plotted as a function of 8 at 5 
different redshifts spanning the overlap epoch (for this model). Lower Panels: The corresponding examples of skewness functions. Two 
cases are shown, corresponding to different values for the critical overdensity prior to the overlap epoch (A c = 5, left panels; and A c = 20, 
right panels). 



points are based on Storrie-Lombardi et al. (1994). Again 
the model is in good agreement with the available observa- 
tions. We note that the observed mean-free-path is found 
from the number density of Ly-limit systems and is inde- 
pendent of the Lya forest absorption derived quantities of 
ionization rate and volume averaged neutral fraction, as well 
as being independent of the HI mass-density measurements. 
Our simple model therefore simultaneously reproduces the 
evolution of three independent measured quantities. 

At a specified redshift, our model yields the mass- 
averaged fraction of ionised regions within the IGM on var- 
ious scales J? as a function of overdensity. We may then cal- 
culate the corresponding 21cm brightness temperature con- 
trast 



3.1 Fluctuations in 21cm emission 

Given the distribution of 8 from the primordial power- 
spectrum of density fluctuations, we may find the proba- 
bility distribution dP/dT of brightness temperature T in 
redshifted 21cm intensity maps, and hence the second and 
third moments of the distribution as functions of angular 
scale and redshift (Wyithe & Morales 2007). The structure 
of our semi-analytic model for reionisation makes it natural 
to discuss the moments of the real space intensity fluctua- 
tions smoothed within top-hat window functions of angular 
radius 6. 

The second moment of these distributions corresponds 
to the auto-correlation function of brightness temperature 
smoothed on an angular radius 6: 



T(S,R) = 22mK 



l + z 
7.5 



1/2 



(l-QiF M (A h S,R)) ( 1 + -8 



(8) 

where the pre-factor of 4/3 on the overdensity refers to the 
spherically averaged enhancement of the brightness tem- 
perature due to peculiar velocities in overdense regions 
(Bharadwaj & Ali 2005; Barkana & Loeb 2005). In the lower 
right panel of Figure Q] we show the evolution of the mean 
brightness temperature contrast of 21cm emission with red- 
shift. For comparison, we also show the fraction of IGM yet 
to overlap. Significant mean redshifted 21cm emission will 
extend for at least a redshift unit beyond overlap (during 
which time it has dropped by a factor of ~ 5). This implies 
that a global step signature will be very gradual, making its 
detection by instruments with a limited band-pass challeng- 
ing. 



<(T-<T» 2 > = 
1 



2na(R) 



d8(T{8,R) - {T)Y e 



where 



(T) 



2ttg{R) 



1 



dS T(8,R)e~ 



(9) 



(10) 



and a{R) is the variance of the density field (at redshift z) 
smoothed on a scale R. 

A measure of the departure of the statistics from those 
of a Gaussian random field is provided by the third moment 
of the distribution dP/ dT. We refer to the third moment as 
the skewness function of brightness temperature smoothed 
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Figure 3. Examples of auto-correlation and skewness functions. Two cases are shown, corresponding to different values for the critical 
overdensity prior to the overlap epoch (A c = 5, left panels; and A c = 20, right panels). Upper Panels: The auto-correlation is plotted as 
a function of redshift for 3 different co-moving scales of relevance for observations using upcoming low-frequency radio telescopes. The 
heavy lines show our model. For comparison the thin lines show the fluctuations in the absence of galaxy bias. Also shown for comparison 
is the fraction of IGM yet to overlap (thick grey line dropping to zero at z = 6 — 7). Lower Panels: The corresponding examples of 
skewness functions. 



with top-hat windows of angular radius 6: 



{(T-{T}Y) = 
1 



2-ko-{R) 



dS (T(S, R) 



<T»V 



(11) 



Examples of auto-correlation and skewness functions 
are plotted in Figure [5] The upper panels show the auto- 
correlation plotted as a function of 6 at 5 different redshifts 
spanning the overlap epoch (for this model). The lower pan- 
els show the corresponding examples of skewness functions. 
Further examples of auto-correlation and skewness functions 
are plotted in Figure [3] but this time as a function of redshift 
for 3 different co-moving scales of relevance for upcoming 
low-frequency radio telescopes. The thick dark lines show 
the evolution for our full model. We also show (thin dark 
lines) the fluctuations that would be present in the absence 
of galaxy bias. Comparison of these curves demonstrates 
that galaxy bias is important in setting fluctuation levels 
prior to the overlap epoch, but becomes progressively less 
important after the overlap epoch is completed. This is be- 
cause the mean-free-path becomes large following the com- 
pletion of overlap, and as a result the ionizing photons that 
form the background at a particular location in the IGM 
are then drawn from a volume that is much larger than the 
regions among which the fluctuations are measured. Galaxy 
bias remains important until lower redshifts on larger scales 
R because the mean-free-path exceeds R only at later times. 
Also shown for comparison (heavy grey lines) is the fraction 
of IGM yet to overlap (1 - Qi). 

Figures [2]|3] suggest that a characteristic feature of the 
overlap epoch will be the inversion of the 21cm intensity dis- 



tribution from negative to positive skewness. As mentioned 
above, the ionization state of an overdense region is set by 
the competing effects of galaxy bias, and higher gas den- 
sity (and recombination rate). With respect to 21cm fluc- 
tuations, this inversion therefore corresponds to the shift in 
dominance from the biased radiation field prior to the over- 
lap epoch, to dominance of gas density and recombinations 
on small scales after the overlap epoch (when the long ion- 
izing mean-free path washes out the effect of galaxy bias). 

To investigate the sensitivity of forthcoming low- 
frequency telescopes to the fluctuations in 21cm emission 
after reionization, we have estimated the signal-to-noise for 
the Mileura-Widefield Array (MWA0) using the formalism 
outlined in Wyithe & Morales (2007). The MWA, which is 
currently under construction will comprise a phased array of 
500 tiles (each tile will contain 16 cross-dipoles) distributed 
over an area with diameter 1.5km. In Figure [4] we show ex- 
amples of signal-to-noise for auto-correlation and skewness 
functions measured with the MWA assuming 1000 hours in- 
tegration (the examples correspond to the models in Fig- 
ure [3]). In the upper panels we plot the signal-to-noise for 
the auto-correlation as a function of redshift for 3 different 
co-moving scales of relevance for the MWA. Also shown for 
comparison as before is the fraction of IGM yet to overlap 
(1 — Qi). In the lower panels we show the corresponding 
signal-to-noise for the skewness function. The model sug- 
gests that the auto-correlation will reach a minimum around 
a redshift unit following the completion of the overlap epoch, 



1 see http://www.haystack.mit.edu/ast/arrays/mwa/index.html 
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Figure 4. Examples of signal-to-noise for auto-correlation and skewness functions measured with the MWA assuming 1000 hours 
integration. Two cases are shown, corresponding to different values for the critical overdensity prior to the overlap epoch (A c = 5, left 
panels; and A c = 20, right panels). Upper Panels: The signal-to- noise for the auto-correlation function is plotted as a function of redshift 
for 3 different co-moving scales of relevance for observations using upcoming low-frequency radio telescopes. Shown for comparison is 
the fraction of IGM yet to overlap according to the definition in § [2] (thick grey lines). Lower Panels: The corresponding examples of 
signal-to-noise for skewness functions. Horizontal lines have been drawn at a signal-to-noise of 3 to guide the eye. 



but that even here it will remain detectable at a signal-to- 
noise of several. Moreover the signal-to-noise at z ~ 3 — 4 
should be comparable with that near the peak of the reion- 
ization era due to the decreased level of sky-noise at lower 
redshifts. The MWA is designed to be sensitive to frequen- 
cies as high as 300MHz, or redshifts as low as z ~ 3.5. Thus 
redshifted 21cm observations could be used to track the en- 
tire overlap era into the post reionization IGM. 



where P(k) is the primordial power-spectrum of the density 
field, extrapolated linearly to z = 0. The bias b^i, which has 
dimensions of temperature squared may be determined di- 
rectly from our calculation of the auto-correlation function, 
i.e. 



621 {Rf 



((T-(T)) 2 ) 



(13) 



4 THE POWER-SPECTRUM OF 21CM 

FLUCTUATIONS AFTER REIONIZATION 

Our analysis thus far has concentrated on the auto- 
correlation function of 21cm fluctuations. A more powerful 
statistical probe will be provided by the power-spectrum 
which is naturally accessible to interferometric observa- 
tions such as those to be carried out by the MWA. Dur- 
ing the reionization epoch, the relation between the power- 
spectrum of 21cm fluctuations and the underlying matter 
power-spectrum is complex, and in the late stages is domi- 
nated by the formation of large ionized bubbles (Furlanetto 
et al. 2004; McQuinn et al. 2006). After reionization (i.e. 
z < 5.5) there are no longer separate ionized bubbles as 
most of the IGM is ionized. Our modeling suggests that the 
skewness of the 21cm intensity fluctuation distribution will 
be small during the post overlap epoch. As a result the 21cm 
emission will follow the density field up to a constant of pro- 
portionality whose value reflects the neutral fraction and a 
bias like factor due to the varying rate of recombination 
with overdensity. Thus, after the overlap of HII regions we 
may write the following expression for the power-spectrum 
of 21cm fluctuations 

ffei(*0 =b 2 21 P(k)D{z) 2 , (12) 



The bias 621 is a function of R, but at z < 5 is nearly 
independent of scale over the range of interest. We show 
examples of the resulting power-spectrum (thick solid lines) 
in Figure[5]at z — 3.5 (left panels) and z — 5 (right panels). 
The model power-spectrum shown was computed assuming 
A c = 20, for which we find 621 ~ 1.2mK at both z — 3.5 
and z — 5. 

Calculations of the sensitivity to the 21cm power- 
spectrum for an interferometer have been presented by a 
number of authors. We follow the procedure outlined by 
McQuinn et al. (2006), drawing on results from Bowman, 
Morales & Hewitt (2006) for the dependence of the array 
antenna density on radius p(r). The uncertainty in a mea- 
surement of the power-spectrum has two separate compo- 
nents, due to the thermal noise of the instrument (AP21,n), 
and due to sample variance within the finite volume of the 
survey (AP2i,sv)- The contamination of foregrounds pro- 
vides an additional source of uncertainty in the estimate 
of the power-spectrum. McQuinn et al. (2006) have shown 
that it should be possible to remove the power due to fore- 
grounds to a level below the cosmological signal, provided 
that the region of band-pass from which the power-spectrum 
is estimated is substantially smaller than the total band- 
pass available. Following the approximation suggested in 
McQuinn et al. (2006), we combine the above components to 
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Figure 5. The power-spectrum of 21cm fluctuations after reionization (thick solid lines) at z = 3.5 (left panels) and z = 5 (right panels). 
The model power-spectrum shown was computed assuming A c = 20. Also shown for comparison are estimates of the noise for MWA 
(upper panels) and its future extension MWA5000 (lower panels). In each case we plot the sample- variance (dotted lines) and thermal 
noise (dashed lines) components of the uncertainty within fc-space bins of size Afc = fc/10. The combined uncertainty including the 
minimum fc cutoff due to foreground subtraction, is also shown as the thin solid lines. 



yield the uncertainty on the estimate of the power-spectrum 

AP 2 i = AP 2 i,sv + AP 2 i,n if k > fc min 

= oo otherwise (14) 

where fc m i n = 2n/y and y is the co- moving line-of-sight 
distance corresponding to the band-pass within which the 
power-spectrum is measured. 

Estimates of the sample-variance (dotted lines) and 
thermal noise (dashed lines) components of the uncertainty 
for detection by the MWA are plotted in the upper panels 
of Figure [5] We model the antennae distribution as hav- 
ing p(r) oc r~ 2 with a maximum radius of 750m and a finite 
density core of radius 20m, and we assume a lOOOhr integra- 
tion, a bandpass of B = 8MHz, and fc-space bins of widtlfl 
Afc = fc/10. The combined uncertainty including the mini- 
mum fc cutoff due to foreground subtraction is shown as the 
thin solid line. The results confirm our expectations based 
on results in the previous section. The power-spectrum of 
21cm fluctuation should be detectable at high significance 
at z ~ 3.5, well after the completion of reionization. In- 
deed the signal-to-noise for detection of the power-spectrum 
at z ~ 3.5 is not dissimilar to expectations for the signal- 
to-noise during the reionization era itself (e.g. McQuinn et 
al. 2006). On the other hand, due to the increasing sky dom- 



2 The signal-to-noise is increased in proportion to V Afc, and so 
will be substantially better per bin in measurements of the power- 
spectrum at lower resolution in k. 



inated system-noise at higher redshift, a detection will be 
more challenging at redshifts just after overlap. 

At values of fc ~ a few xlO" 1 Mpc" 1 , the measure- 
ment of the power-spectrum using the MWA will be limited 
by the thermal sensitivity of the array, and so the signal- 
to-noise achievable in this regime will be greatly enhanced 
by a subsequent generation of instruments with larger col- 
lecting area. As an example we consider a hypothetical fol- 
lowup instrument to the MWA which would comprise 10 
times the total collecting area. We refer to this followup 
telescope as the MWA5000. The design philosophy for the 
MWA5000 would be similar to the MWA, and we therefore 
assume antennae distributed as p(r) oc r~ 2 with a diameter 
of 2km and a flat density core of radius 80m (see McQuinn 
et al. 2006). In the lower panels of Figure [5] we present esti- 
mates for measurement of the 21cm power-spectrum using 
MWA5000. The sample-variance (dotted lines) and thermal 
noise (dashed lines) components of the uncertainty are plot- 
ted as before assuming a lOOOhr integration, a bandpass of 
B = 8MHz, and fc-space bins of width Afc = fc/10. The 
combined uncertainty including the minimum fc cutoff due 
to foreground subtraction is shown as the thin solid line. 
This figure, in combination with previous calculations of 
sensitivity to the 21cm power-spectrum during reionization 
(e.g. McQuinn et al. 2006) demonstrates that the MWA5000 
could measure the 21cm power-spectrum in detail across the 
entire frequency range accessible to an MWA antennae, cor- 
responding to a redshift range of 3.5 < z < 15. 

The reionization epoch has formed the dominant driver 
for studies of the statistical properties of redshifted 21cm 
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emission to date, while we have argued in this paper that 
the study of the IGM using 21cm emission could also be 
extended to lower post-reionization redshifts. However the 
measurement of a 21cm power spectrum will also provide 
significant cosmological information (McQuinn et al. 2006). 
For example the angular scale of the baryonic acoustic peak 
in the matter power spectrum, measured through redshifted 
21cm observations could be used to constrain the evolu- 
tion of the dark energy with cosmic time (Blake & Glaze- 
brook 2003; Seo & Eisenstein 2007; Angulo et al. 2007) in 
the redshift range of 3.5 < z < 6, as well as during the 
reionization era. No other probes are currently effective dur- 
ing these epochs (Corasaniti, Huterer, & Melchiorri 2007). 
A detailed analysis of the prospects for such a study with 
future low-frequency arrays will be presented elsewhere. 



5 DISCUSSION 

While the recent discovery of distant quasars has allowed 
detailed absorption studies of the state of the IGM at z ~ 6 
(Fan et al. 2006; White et al. 2003), the interpretation of 
the data is not yet robust. The prime reason for the ambi- 
guity in interpreting high redshift quasar absorption spec- 
tra is that Lya absorption can only be used to probe vol- 
ume averaged neutral fractions that are smaller than 10 -3 , 
owing to the large cross-section of the Lya resonance. A 
better probe of the process of reionization is provided by 
redshifted 21cm observations. The conventional wisdom is 
that the 21cm signal disappears after the overlap epoch, be- 
cause there is little neutral hydrogen left through most of 
the intergalactic space. But observations of damped Lya sys- 
tems out to redshift z ~ 4 show the density parameter of HI 
to be JIhi ~ 10 -3 , indicating that the mass averaged neu- 
tral hydrogen fraction remains at the level of a few percent 
throughout cosmic history. 

Existing constraints from Lya studies are based on line- 
of-sight averaged optical depths and are therefore sensitive 
to the volume filling fraction of neutral hydrogen. These 
observations are limited in their ability to probe the IGM 
once the overlap era is approached due to the high opti- 
cal depth to Lya absorption. On the other hand, redshifted 
21cm observations directly probe the pockets where most of 
the hydrogen resides after reionization. In this paper we have 
considered the detection of fluctuations in 21cm emission 
from neutral hydrogen after the completion of reionization 
using low-frequency telescopes that are currently under con- 
struction (specifically the Mileura Widefield Array, MWA). 
Based on a physically motivated model for the reionization 
era we have shown that the residual neutral gas, which is 
contained in high density pockets following the end of the 
overlap epoch, produces a substantial signal of 21cm fluctua- 
tions. These fluctuations directly probe the post-reionization 
IGM and have a power-spectrum that will be detectable by 
the MWA at a signal-to-noise ratio that is comparable to 
observations during the reionization epoch. We stress that 
in the case of observations at z < 5 the prediction of the 
power-spectrum amplitude is quite robust, being primar- 
ily dependent on the underlying mass power-spectrum and 
the mass-averaged neutral fraction of hydrogen which, un- 
like predictions for the reionization epoch, represent well- 
measured quantities. 



Models of reionization predict that non-Gaussianity in 
the distribution of 21cm intensities would be the signature 
of star-formation driven reionization of the IGM. Our cal- 
culations predict that after the overlap of HII regions in the 
IGM, there will be a transition between an epoch of nega- 
tive skewness when the statistics of the intensity distribution 
are driven by galaxy- bias (prior to overlap) , and an epoch of 
positive skewness when they are driven by the density distri- 
bution and recombination rate (post-overlap). Observation 
of this rapid transition in the statistics of the 21cm signal 
will therefore provide added confidence that the end of the 
reionization era has been detected. 

Observations of the post-reionization 21cm power- 
spectrum will be significantly enhanced by 2nd generation 
low-frequency telescopes with increased collecting area. In 
addition to the utility of 21cm fluctuations for studies of 
the evolution in the ionization state of the IGM during the 
post reionization epoch, measurement of the 21cm power- 
spectrum may therefore yield precise measurements of the 
angular scale of the acoustic peak in the matter power- 
spectrum. Thus in the future it may become possible to use 
redshifted 21cm studies to constrain the dark energy in the 
unexplored redshift range of 3.5 < z < 15, where it is not 
accessible through other techniques. 
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